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Abstract 

Thermodynamical Behaviour of Inhomogeneous Universe with Varying A 
in Presence of Electromagnetic Field is obtained. F\2 is the non-vanishing 
component of electromagnetic field tensor. To get a deterministic so- 
lution, it is assumed that the free gravitational field is Petrov type-II 
non-degenerate. The value of cosmological constant is found to be small 
and pasitive supported by recent results from the supernovae observations 
recently obtained by High-Z Supernovae la Team and Supernovae Cosmo- 
logical Project. A relation between cosmological constant and thermody- 
namical quantities is established. Some physical and geometric properties 
of the model are also discussed. 
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1 Introduction 

A lesson given by the history of cosmology is that the concept of the cosmolog- 
ical term revives in the days of crisis and we have more reasons than ever to 
belive that the cosmological term is the necessary ingredient of any cosmological 
model. The possibility of adding a cosmological vacuum energy density to the 
Einstein field equations raises the question empirical justification of such a step. 
A positive cosmological constant helps overcome the age problem, conected on 
the one side with the high estimates of the Hubble parameters and with the age 
of the globular clusters on the other. Further, it seems that in order to retain 
the cold dark matter theory in the spacially flat universe most of the critical 
density should be provided by a passitive cosmological constant [IJ [2] . Observa- 
tionnal data indicate that the cosmologial constant, if nonzero, is smaller than 
10~ 55 cm~2. However, since everything that contributes to the vacuum energy 
acts as a cosmological constant it can not just be dropped without serious con- 
siderations. Moreover particle physics expectations for A exceeds its present 
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value by the factor of order 10 120 ie in a sharp contrast to observations. To ex- 
plain this apparent discrepancy the point of view has been adopted which allow 
the A term to vary in time [3] — [TT] . The idea of that during the evolution of 
universe the energy density of the vacuum decays into the particles thus leading 
to the decrease of the cosmological constant. As the result one has the creation 
of particles although the typical rate of the creation is very small. 

The creation of matter and entropy from vacuum has been studied via quan- 
tum field theory in curved spacetime [T^l [T3] . Most cosmological models exhibit 
a singurarity which presents difficulties for interpreting quantum effects, because 
all macroscopic parameters of created particles are infinite there. This leads to 
the problem of the initial vacuum. A regular vacuum for species of the created 
particles can be defined in simple terms as a state where all mean values de- 
scribing the particles, such as energy density, number density, entropy etc are 
zero. But this simple condition is not achieved in many scenarios, so that cither 
one has to postulate an initial state beyond the singularity, or to assume that 
there was a non-zero number of particles at the initial vacuum. One attempt 
to overcome these problems is via incorporating the effect of particle creation 
into Einstein's field equactions. In the present study I interpret the source 
of created particles as a decaying vacuum, described phcnomenologically by a 
time-dependent cosmological constant A(i). 

There are significant observational evidence for the detection of Einsteins 
cosmological constant, A or a component of material content of the universe that 
varies slowly with time and space to act like A. Some of the recent discussions 
on the cosmological constant problem and on cosmology with a time-varying 
cosmological constant by Ratra and Peebles [14], and Sahni and Starobinsky 
[15] , point out that in the absence of any interaction with matter or radiation, 
the cosmological constant remains a constant. However, in the presence of in- 
teractions with matter or radiation, a solution of Einstein equations and the 
assumed equation of covariant conservation of stress-energy with a time- varying 
A can be found. This entails that energy has to be conserved by a decrease in the 
energy density of the vacuum component followed by a corresponding increase 
in the energy density of matter or radiation (see also Carroll, Press and Turner 
[15] , Peebles [T7J, Padmanabhan [T5]). There is a plethora of astrophysical ev- 
idence today, from supcrnovac measurements (Pcrlmutter et al. 19J, Riess et al. 
[20] . Garnavich et al. [21] . Schmidt et al. [22], Blakeslee et al. [23], Astier et al. 
[24]), the spectrum of fluctuations in the Cosmic Microwave Background (CMB) 
[25] . baryon oscillations [26] and other astrophysical data, indicating that the 
expansion of the universe is currently accelerating. The energy budget of the 
universe seems to be dominated at the present epoch by a mysterious dark en- 
ergy component, but the precise nature of this energy is still unknown. Many 
theoretical models provide possible explanations for the dark energy, ranging 
from a cosmological term [37] to super-horizon perturbations [33] and time- 
varying quintessence scenarios |29| . These recent observations strongly favour 
a significant and a positive value of A with magnitude A (^r) ~ 10 -123 . 

The standard Friedman-Robertson- Walker (FRW) cosmological model pre- 
scribes a homogeneous and an isotropic distribution for its matter in the de- 
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scription of the present state of the universe. At the present state of evolution, 
the universe is spherically symmetric and the matter distribution in the universe 
is on the whole isotropic and homogeneous. But in early stages of evolution, 
it could have not had such a smoothed picture. Close to the big bang singu- 
larity, neither the assumption of spherical symmetry nor that of isotropy can 
be strictly valid. So we consider plane-symmetric, which is less restrictive than 
spherical symmetry and can provide an avenue to study inhomogeneities. Inho- 
mogeneous cosmological models play an important role in understanding some 
essential features of the universe such as the formation of galaxies during the 
early stages of evolution and process of homogenization. The early attempts 
at the construction of such models have done by Tolman [3U] and Bondi [3"T] 
who considered spherically symmetric models. Inhomogcneous plane-symmetric 
models were considered by Taub [32[ [33] and later by Tomimura [34] , Szekeres 
[35] . Recently, Senovilla [36^ obtained a new class of exact solutions of Einstein's 
equation without big bang singularity, representing a cylindrically symmetric, 
inhomogeneous cosmological model filled with perfect fluid which is smooth and 
regular everywhere satisfying energy and causality conditions. Later, Ruis and 
Senovilla [37] have separated out a fairly large class of singularity free mod- 
els through a comprehensive study of general cylindrically symmetric metric 
with separable function of r and t as metric coefficients. Dadhich et al. [35] 
have established a link between the FRW model and the singularity free fam- 
ily by deducing the latter through a natural and simple in-homogenization and 
anisotropization of the former. Recently Bali and Tyagi[39 , Pradhan et al 40J 
obtained a plane-symmetric inhomogeneous cosmological models of perfect fluid 
distribution with electro-magnetic field. 

The occurrence of magnetic fields on galactic scale is well-established fact 
today, and their importance for a variety of astrophysical phenomena is generally 
acknowledged as pointed out by Zeldovich et al. 41J. Also Harrison [42] has 
suggested that magnetic field could have a cosmological origin. As a natural 
consequences, we should include magnetic fields in the energy-momentum tensor 
of the early universe. The choice of anisotropic cosmological models in Einstein 
system of field equations leads to the cosmological models more general than 
Robertson- Walker model [13]. Strong magnetic fields can be created due to 
adiabatic compression in clusters of galaxies. Primordial asymmetry of particle 
(say electron) over antiparticle (say positron) have been well established as C 
P (charged parity) violation. Asseo and Sol [44] speculated the large-scale inter 
galactic magnetic field and is of primordial origin at present measure 10~ 8 G 
and gives rise to a density of order 10 -35 gcm~ 3 . The present day magnitude of 
magnetic energy is very small in comparison with the estimated matter density, 
it might not have been negligible during early stage of evolution of the universe. 
FRW models are approximately valid as present day magnetic field is very small. 
The existence of a primordial magnetic field is limited to Bianchi Types I, II, 
III, VIq and VIIq as shown by Hughston and Jacobs [IS]. Large-scale magnetic 
fields give rise to anisotropies in the universe. The anisotropic pressure created 
by the magnetic fields dominates the evolution of the shear anisotropy and it 
decays slower than if the pressure was isotropic [46] [47]. Such fields can be 
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generated at the end of an inflationary epoch [H] — [SU]. Anisotropic magnetic 
field models have significant contribution in the evolution of galaxies and stellar 
objects. 



2 The metric and field equations 

We consider the metric in the form of Marder [51] 

ds 2 = A 2 {dx 2 - dt 2 ) + B 2 dy 2 + C 2 dz 2 , 



(1) 



where the metric potential A, B and C are functions of x and t. The energy 
momentum tensor is taken as 



T i — (P + P) v i vj + P9i + E l ! 

where Ef is the electro-magnetic field given by Lichnerowicz 



as 
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hih'fav* + ^g 3 i )-h i h j 



(2) 



(3) 



Here p and p are the energy density and isotropic pressure respectively and v 1 
is the flow vector satisfying the relation 



ffy«V = -1. 

p, is the magnetic permeability and hi the magnetic flux vector defined by 

hi — FjiV^ , 
A* 

where *Fy is the dual electro-magnetic field tensor defined by Synge [55] 



" F- 



(4) 



(5) 



(6) 



Fij is the electro-magnetic field tensor and tijki is the Levi-Civita tensor density. 
The coordinates arc considered to be comoving so that v 1 = = v 2 = v 3 and 
v 4 = -j. We consider that the current is flowing along the z-axis so that ^ 0, 
h\ = = hi = hi- The only non- vanishing component of Fij is Fx%. The 
Maxwell's equations 

Fij. k + F jk ;i + Fki-j = 0, (7) 



and 



0. 



(8) 



require that F\2 be function of x alone. We assume that the magnetic perme- 
ability as a function of x and t both. Here the semicolon represents a covariant 
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differentiation. 



The Einstein's field equations ( in gravitational units c = 1, G = 1 ) read as 
l4-±Rgi+Agi = -8wTi, (9) 
for the line element (1) has been set up as 
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^TT + ^-Tllf + ^J-Tlt + ; ?J' < 14 » 

where the sub indices f and 4 in A, B, C and elsewhere indicate ordinary dif- 
ferentiation with respect to x and t, respectively. 



3 Solution of the field equations 

Equations (JTUJ) - dHJ) lead to 

A4 \ ^.4 / C4 \ B^C^ 



Ai 
and 



A B A \ B C ) BC 

a (constant), (15) 

^ (16) 

Eqs. (HHJ) - HU) represent a system of five equations in seven unknowns A, 
B, C, p, p, A and p. For the complete determination of these unknowns two 
more conditions are needed. As in the case of general-relativistic cosmologies, 
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the introduction of inhomogeneities into the cosmological equations produces a 
considerable increase in mathematical difficulty: non-linear partial differential 
equations must now be solved. In practice, this means that we must proceed 
either by means of approximations which render the non-linearities tractable, 
or we must introduce particular symmetries into the metric of the space-time 
in order to reduce the number of degrees of freedom which the inhomogeneities 
can exploit. In the present case, we assume that the metric is Petrov type-II 
non-degenerate. This requires that 
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Let us consider that 



A = f(x)X(t), 
B = g(x)n(t), 
C = g{x)v{t). 
Using (JUS) in (pl|) and CLZ|), we get 
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where m and n are constants of integration. Equations (fT5|) . IT8|) and ((20|) lead 
to 

U4I/4 



Equation (|19p leads to 
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Equations fl20J), ([2SJ and (QU) lead to 

W± = Me Lt+2 ^ u , 
where M is constant. From equations [[25]). 1(25]). (gBJ) and (|27|) , 

(6 - l)C/ 44 + 2(6 - l)C/f - 26Afe Lt+2 ( b - 1 ) c/ C/ 4 - MLe Lt+ ^ b ~ 
If we put e 2U = £ in equation (|28|) . we obtain 



If we consider £ = e ht , then equation (|29|) leads to 

(*zl) ff 2 e « _ Af-(e Lt e' ifct ) = ae' 14 , 
which again reduces to 

h L 



and 
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Equations (|2"T)) and reduce to 

iy = Aft + log TV, 
where N is an integrating constant. Eq. (|2^|) leads to 

5 = /Jsinh^TT (ai + 5), 

where 



= Po ( 



(2-6) ' 

and /3o, S being constants of integration. Hence 

b-2 

/ = n/3sinh 2 (''- 1 ) (aar + 5), 
A = me^^ , 



7 



Therefore, we have 

A = fX = mnPe*ft=V sinh^ (ax + S), (40) 
B = gfj, = jV/3e(^ +2M H sinh 3 ^) (ax + 5), (41) 
C = gv= -^e(T^" 2M )* sinh^T (ax + 5). (42) 
By using the transformation 

aX = ax + 8, 

Y = Gy, 
Z = Hz, 

T = t, (43) 

the metric (JT]) reduces to the form 

ds 2 = K 2 sinh b - 2 (aX)e^ (dX 2 - dT 2 )+ 

sinh^ (aX)e^ +2M ) T dY 2 + sinh^ (aX)e^~ 2M ) T dZ 2 , (44) 
where K = mn(3, G = N/3 and H = 

4 Some Physical and Geometric Features 

The physical parameters, pressure (p) and density (p), for the model (|44|) are 
given by 



8np = —2 e^r sinn 2 " 6 (qX) 
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In this case to find the explicit value of cosmological constant A(i), one may 
assume that the fluid obey an equation of state of the form 



p = 7/5 



(47) 
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Figure 1: The plot of energy density p vs. time T with parameters 6 — 0.5, m = 
0.03, n = 0.25, K = 0.2 and 7 = 0.4 



where 7(0 < 7 < 1) is a constant. 

Using equation (|4"T)l in (|4"5)) and then solving with (|4^|) . we have 
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Eliminating 7 from (|4l))) and l|48p. we obtain 
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From (|4^|) and (T4"8"|) , we note that p(t) is the decreasing function of time and 
p > for all time. This behaviour is clearly shown in figure 1, as a representa- 
tive case with appropriate choice of constants of integration and other physical 
parameters using reasonably well known situations. 
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Figure 2: The plot of cosmological constant A vs. time T with parameters 
b = 0.5, m = 0.03, n = 0.25, K = 0.2 and 7 = 0.4 



In spite of homogeniety at large scale our universe is inhomogenrous at small 
scale, so physical quantities being position dependent are more natural in our 
observable universe if we do not go to super high scale. This result show this kind 
of physical importance. In recent time theoreticians and observers draw high 
attention on A - terms due to various reasons. The non-trivial role of the vacuum 
in early universe generates a A that leads to inflationary phase. Obsevationally 
this term provides an additional parameter to accommodate conflicting data on 
the value of Hubble's constant, deceleration parameter, the density parameters 
and the age of universe [54j and |55j . Assuming that A owes it origin to vacuum 
interactions, as suggested in particular by Sakharav [56 , it follows that it would 
in general be a function of space and time co-ordinates, rather than a strict 
constant. In a homogeneous universe A will be most time dependent |57j . In 
our case this approach can generate A that varies both with space and time. In 
considering the nature of local massive objects, however the space dependence 
of A can not ignored. For detail discussion, the readers are advised to see the 
reference (Narlikar, Pecker and Vigier [55], Ray and Ray [55], Tiwari, Ray and 
Bhadra [HO]). 

From ([49]) . we see that cosmological constant A is a decreasing function of 
time and it approaches a small positive value at late time. This behaviour is 
clearly shown in figure 2. Recent cosmological obsevations suggest the existence 
of positive cosmological constant A with the mag nitudeA(^) » lO" 123 . These 
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observations on magnitude and red-shift of type la supernova suggest that our 
universe may be an accelerating one with induced cosmological density through 
the cosmological A-term. Thus the model presented in this paper is consistent 
with the results of recent observations. 

The non-vanishing component F12 of the electromagnetic field tensor is given 
by 



= t lfi2ML (A+ 
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where [i remains undetermined as function of x and t both. 

The scalar of expansion {&) calculated for the flow vector (V) is given by 
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The shear scalar (a 2 ), acceleration vector (ii), deceleration parameter q, proper 
volume V and Hubble parameter H for the model (|4"5)) are given by 
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From equations ([51]) and (f52"| . we have 
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The rotation u is identically zero and the non- vanishing component of conformal 
curvature tensor are given by 
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The reality conditions (Ellis [BT] ) 

(i)p + p>0, (ii)p + 3p> 0, 

lead to 
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respectively. 

The dominant energy condition is given by Hawking and Ellis [62 
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5 Thermodynamical behaviour and entropy of 
universe 

From the thermodynamics [63 , 64 , we apply the combination of first and second 
law of thermodynamics to the system with volume V. As we know that 

TdS = d{pV) +pdV (66) 
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where T, S represents the tempreture and entropy respectively. 
Eq. (|66|) may be written as 

TdS = d[(p + p)V] - Vdp (67) 

The integrability condition is necessary to define a perefect fluid as a thermo- 
dynamical syetem |65l 166] . It is given by 

dp = ^^rfT (68) 

Plugging eq. (|68|) in eq. (167|) . we have the differential equation 

dS=^d[(p + p)V]-(p+p)V^ (69) 



we rewrite eq (fBU)) as 

'(p + p)V 
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where c is constant. 

Hence the entropy is defined as 



T 

Let the entropy density be s, so that 



S=?±IV (71) 
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where p — jp and < 7 < 1. 

If we define the entropy density in terms of temprature then the first law of 
thermodynamics may be written as 

d(pV)+ 1P dV=(l + p)Td(^j (73) 

which on integration yields 

T = CQpT^TT) (74) 
where Co is constant of integration. 

From eqs. (|72|) and (71), we obtain 

'1 + 7 



Co 



(75) 



These equation are not valid for 7 = — 1. For the Zel'dovich fluid (7 = 1), we 
get 

T = c pi (76) 
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Thus the entropy density is proportional to the tempreture. 

Now the tempreture, entropy density and entropy of Zcl'dovich universe is given 

by 



T = T n 
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where T = s 



— %= and 5*0 = snK are constant. 
For radiating fluid (7 = |), we get 
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Thus the entropy density is proportional to cube of tempreture. 

Now the tempreture, entropy density and entropy of radiating universe is given 

by 
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where T 00 = (87r c °i /2 , sqo = E^pfpTa an( i -Sbo = s K 2 are constant. 



6 Discussion and Concluding Remarks 

The present study deals the plane-symmetric inhomogeneous cosmological model 
of electro-magnetic perfect fluid as the source of matter. FRW models are ap- 
poximatly valid, as present day magnetic field is very small. Maarteens [67] in 
his study explained that magnetic fields are observed not only in stars but also in 
galaxies. In princple, these fields could play a significant role in structure forma- 
tion but also affect the anisotropies in cosmic microwave background radiation 
[CMB]. Since the electric and magnetic fields are interrelated, their indepen- 
dent nature disappears when we consider them as time dependance. Hence, it 
would be proper to look upon these fields as a single field - electromagnetic field. 
Generally the model represents expanding, shearing, non-rotating and Pctrov 
type-II non-degenerate universe in which the flow vector is geodesic. We find 
that the model starts expanding at T = and goes on expanding indefinitely. 
For large values of T, the model is conformally flat and Petrov type-II non- 
degenerate. Since f = constant, hence the model does not approach isotropy. 
The value of cosmological constant is found to be small and pasitive supported 
by recent results from the supernovae observations recently obtained by High-Z 
Supernovae la Team and Supernovae Cosmological Project. The relation between 
cosmological constant and thermodynamical quantities is established and from 
eqs. (|82|) and (85), it is clear that cosmological constant affect entropy. The 
electromagnetic field tensor does not vanish when L ^ 0, M ^ 0, and 6^1. 
For large values of T and L + 2M(1 — 6) < 0, F%2 tends to zero. For 6 = 1, we 
obtain singularity and model approach isotropy. 6 < 1 and 6 > 1 imposed the 
restriction on the value of M and L which affect all the physical and kinematical 
parameters of the model. 

In spite of homogeneity at large scale our universe is inhomogeneous at small 
scale, so physical quantities being position-dependent are more natural in our 
observable universe if we do not go to super high scale. Our derived model shows 
this kind of physical importance. The expressions for deceleration parameter (q) 
and Hubble parameter (H) given by Eqs. ([54]) and (f56|) respectively are func- 
tions of x and t as in the case of inhomogeneous cosmological models. Obviously 
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these two expressions are different from that of FRW cosmological model. The 
present study also extend the work of Yadav and Bagora [68 for inhomogc- 
neous universe with varying A and clarify thermodynamics of inhomogeneous 
universe by introducing the integrability condition and tempreture. It is found 
that s ~ pa r*j T and s ~ pi ~ T 3 , for Zeldo'vich and radiating fluid model re- 
spectively. From these expressions it is clear that p is the function of tempreture 
and volume, thus a new general equation of state describing the Zel'dovich fluid 
and radiating fluid models as a function of tempreture and volume is found. 
The basic equations of thermodynamics for inhomogeneous universe has been 
deduced which may be useful for better understanding of evolution of universe. 
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